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Abstract. Recent studies have found that the spectral evolution of pulses within 
gamma-ray bursts (GRBs) is consistent with simple radiative cooling. Perhaps more 
interesting was a report that some bursts may have a single cooling rate for the multiple 
pulses that occur within it. We determine the probability that the observed "cooling 
rate invariance" is purely coincidental by sampling values from the observed distribu- 
tion of cooling rates. We find a 0.1 — 26% probability that we would randomly observe 
a similar degree of invariance based on a variety of pulse selection methods and pulse 
comparison statistics. This probability is sufficiently high to warrant skepticism of any 
intrinsic invariance in the cooling rate. 



Much progress has been made in the past few years regarding the spectral evo- 
lution of gamma-ray bursts. Early studies described the time evolution of a single 
hardness parameter in typically fewer than 20 bursts [1-3] . Use of the BATSE LAD 
data has allowed time-resolved spectroscopy of approximately one hundred bursts 
[4]. It has become clear that the time- integrated spectra is not at all representative 
of the time-resolved spectra [5,6]. Several new trends have also been found in anal- 
yses of the time-resolved spectra. In particular, Liang & Kargatis [7] discovered 
a unique relationship between E pea k (the energy at which the vF v spectrum is a 
maximum) and the photon fluence <3>(t) = j t ,l o Fu^df, where 



and $o represents an effective cooling constant. It was also reported that $o ap- 
peared to be invariant from pulse to pulse within some multipulse bursts [7] . Figure 
1 shows an example of this observed invariance in cooling rate for GRB 921207 as 
seen by BATSE LAD 0. 

If the decay of E peak with respect to $ is assumed to be the radiative cooling of 
a plasma, then energy conservation dictates a proportionality between $ an d the 
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FIGURE 1. Evolution of the spectral energy break (squares) and photon flux (histogram) as 
a function of integrated photon flux $ for GRB 921207 as seen by the BATSE LAD 0. The 
cooling decay constant $o is shown for the three observed pulses. The second and third pulses 
are examples of the observed "invariance" of cooling rate seen within gamma-ray bursts. 

number of cooling leptons in the plasma [7]. This in turn suggests that multiple 
pulses observed within a single burst originate from the same plasma. Verification of 
the observed "cooling rate invariance" would have serious implications on gamma- 
ray burst models. A scenario in which a single plasma emits several pulses of 
radiation is difficult to reconcile with many catastrophic burster models [8]. As a 
strict test of the conclusions of Liang & Kargatis [7] we determine the probability 
that the degree of invariance that they observed in their bursts would originate 
from bursts with no intrinsic invariance at all. As described below, this probability 
is high enough to raise doubts about any intrinsic invariance. 



I TESTING THE INVARIANCE OF COOLING RATES 

WITHIN BURSTS 

Liang & Kargatis [7] fitted exponential decays of E peak with respect to fluence 
$ in multipulse bursts and reported that the decay constant <3>o did not change 
from pulse to pulse in many of their bursts. We tested this conclusion by randomly 
sampling values from the observed distribution of $o- Examining the first pulse of 
57 bursts, we found that we can approximate this distribution as $o = 10 L53±a39 
ergs cm -2 . We also assigned a la-confidence region for each random value of $o 
based on the observed distribution cr$ /$ = I0~ a76±0 - 69 . 

In selecting the E pea k decays manually, it is likely that some systematic bias 



TABLE 1. The probabilities of getting such good values for each of 
our three goodness-of-fit statistics G, where the "Best" and "Worst" 
case selections are discussed in the test. Any reasonable analysis 
should fall somewhere between these two cases. These probabilities 
are high enough (> 0.001) to suggest that any invariance seen in the 
data is purely coincidental. 
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occurs favoring the trends we are looking for, namely well fit decays with similar 
decay rates. To account for this, we used two different methods for selecting pulse 
decays. For the case we designate "Best", we selected the decay phases of pulses 
by directly examining plots of E peak versus $. For the case we designate "Worst", 
we instead selected decay phases by examining plots of photon flux versus time. 
In both cases, however, we fit out exponential decay to E peak versus Thus, in 
the "Worst" case scenario, we have blindly selected what portions of the E peak — $ 
plot will be fit. Occasionally, there are some pulses where E peak may actually be 
rising with respect to $ in what was chosen to be a decay phase. We discard these 
instances of "negative decay" . What is left serves as an extreme limit of what our 
final results would have been without our systematic bias. 
We compared two of each bursts' M pulses at a time using 

v2 _ [$o(i)-$ (j)] 2 
ij rr 2 4- a 2 

and then distilled the comparisons within each burst into a single statistic to rep- 
resent that burst (See Table 1 for definitions). Each of these statistics are tailored 
for different null hypotheses. The statistic G\ tests if at least two pulses in a burst 
are similar (and thus "invariant"), while G 2 tests if all the pulses have a similar 
decay constant. G 3 tests for either a single good pairing or several moderately close 
pairings. We believe that this last statistic is the most reasonable for testing our 
results since it does not require that all pulses decay at the same rate (as G 2 does) 
but also does not discard information about multiple pulses repeating (as G\ does). 
Finally, we calculated a table of probabilities ~P(G, M) for our goodness-of-fit statis- 
tics G based on random sampling from our distributions of $ and a <s> /^o- We 
define P(G, M) to be the probability of randomly getting a value of G or lower if we 
compare M pulses drawn randomly from the global distribution of $ - Examining 
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the distibution of P values for our sample of 40 "Best" multipulse bursts, we found 
K-S probabilities of 0.1 —6.5% that the observed "repetitions" occur by chance [9]. 
This number is much higher (12 — 26%) if we remove our selection systematics. 
These high probabilities suggest that the observed invariance in cooling rate is not 
physical. 

II DISCUSSION 

An invariant cooling rate for pulses in multipulse gamma-ray bursts is certainly 
an exciting possibility. Proof of such a pattern would place hard limits on burst 
emission mechanisms. Liang and Kargatis [7] found that 7 out of the 12 multipulse 
bursts in their sample were consistent with an invariant decay (difference < la). 
However, using our "Best" selection and the G\ goodness-of-fit measure, we expect 
to randomly get 5.7 ± 1.3 bursts out of 12 consistent with having an invariant 
decay. We conclude that no intrinsic invariance of $ is required to explain the 
observations reported in Liang and Kargatis [7] . The observed "repeating" cooling 
rates are consistent with being randomly selected from the observed distribution 
of $ values. While these new results do not necessarily rule out cooling rate 
invariance in some bursts, they do show that the observations as an ensemble do 
not require such an invariance. 
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